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1. Introduction 



The existence of Dark Matter seems to require physics beyond the Standard Model. If this physics 
arises from a string/M theory vacuum, one is faced with various problems associated with the 
moduli fields, which are gauge-singlet scalar fields that arise when compactifying string/M theory 
to four dimensions. In particular, moduli fields can give rise to disastrous cosmological effects. 

For example, the moduli have to be stabilized, or made massive, in accord with cosmological 
observations. Even if these moduli are made massive, there could be a large amount of energy stored 
in them leading to the formation of scalar condensates. In most cases, this condensate will scale like 
ordinary matter and will quickly come to dominate the energy density. The moduli are unstable 
to decays to photons, and when this occurs, the resulting entropy can often spoil the successes of 
big-bang nucleosynthesis (BBN). This is the cosmological moduli problem [1-5]. In supersymmetric 
extensions of the standard model, the overproduction of gravitinos can cause similar problems and 
have been a source of much investigation [6-19]. 

In addition, the "standard" picture in which Dark Matter (DM) particles are produced during 
a phase of thermal equilibrium can be significantly altered in the presence of moduli. The moduli, 
which scale like non-relativistic matter, typically dominate the energy density of the Universe 
making it matter dominated. Therefore, the dominant mechanism for production of DM particles is 
non-thermal production via the direct decay of moduli^ . However, this can lead to further problems 
since it is easy to produce too much dark matter compared with what we observe today. 

In recent years there has been considerable progress in our understanding of moduli dynamics 
and their potential in different frameworks which arise in various limits of string/Af theory. The 
most popular examples include the KKLT and Large Volume frameworks in Type IIB string theory 
[23-25], where all moduli are stabilized by a combination of fluxes and quantum corrections. These 
frameworks are also attractive in the sense that they provide a mechanism for supersymmetry 
breaking at low scales TcV), thus accommodating the hierarchy between the Electroweak and 
Planck scales (see [26-28] for reviews). Since one can concretely study the couplings between 
moduli and matter fields, we have an opportunity to address many issues in particle physics and 
cosmology from an underlying microscopic viewpoint. The cosmological moduli/gravitino problems 
and adequate generation of dark matter within the Type IIB frameworks has met with some mixed 
success in a recent paper [29] . 

In this paper we will study a different framework, in which we will also address the Dark Matter 
and modiili/gravitino problems. This is the low energy limit of M theory vacua in which the extra 
dimensions form a manifold of G2 holonomy. Although the study of such vacua has proven to 
be technically challenging, much progress has been made towards understanding the effective four 
dimensional physics emerging from them [31-34]. This includes many phenomenological implica- 
tions of these vacua, in particular relating to issues such as constructing a realistic visible sector 
with chiral matter and non-abelian gauge bosons, supersymmetry breaking, moduli stabilization in 
a dS vacuum as well as explaining the Hierarchy between the Electroweak and Planck scales, as 
exemplified in a number of works [30,35-39]. 

We will show that the moduli, gravitino and dark matter problems are all naturally solved 
within this framework. Because of the presence of moduli, the Universe is matter-dominated from 
the end of inflation to the beginning of BBN. The LSPs are mostly produced non-thermally via 
moduli decays. The final result for the relic density only depends on the masses and couplings of 
the lightest of the moduli (which decay last) and the mass of the LSP. This is related to the fact 
that the LSP is a Wino in the G2-MSSM and that there is a fairly model independent critical LSP 
density at freeze out. For natural/reasonable choices of microscopic parameters defining the G2 

^For other phenomenologically based approaches to non-thermal dark matter and the related issue of baryon 
asymmetry in the presence of scalar decay see [20-22]. 
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framework, one finds that it is possible to obtain a relic density of the right order of magnitude 
(up to factors of 0(1)). With a more sophisticated understanding of the microscopic theory, one 
might obtain a more precise result. The qualitative features which are crucial in solving the above 
problems may also be present in other realistic string/M theory frameworks. 

Moduli which decay into Wino LSPs have been considered previously in the context of Anomaly 
Mediated Supersymmetry Breaking Models (AMSB) by Moroi and Randall [40]. The moduli and 
gaugino masses they consider are qualitatively similar to those of the G2-MSSM. There are some 
important differences however. In particular, the MSSM scalar masses in the [40] arc miich lower 
than the G2-MSSM, leading to much fewer LSPs produced per modulus decay compared to the G2 
models. Furthermore, unlike in AMSB, in the G2 case one is able to calculate all the moduli masses 
and couplings explicitly which leads to a more detailed understanding. In essence, though, many of 
the important ideas in our work arc already present in [40]. The G2-MSSM models can be thought 
of as a concrete microscopic realization of the relevant qualitative features of the AMSB models. 

Interestingly, our actual result for the relic density (equation 6.17) is a few times larger than 
the WMAP value if we use central values for the microscopic constants, which should probably be 
regarded as a success. It is also worth remarking that, contrary to common views, it is not at all 
possible to get any value one wants - we can barely accommodate the actual observed value in the 
G2 framework. 

The paper is organized as follows. In Section 2 we briefly summarize early universe cosmology 
in the presence of moduli, and address many of the issues associated with their stabilization and 
decay. In Section 3 we give a non-technical overview of the main results. This is largely because 
much of this paper involves technical calculations. In section 4 we present a brief review of the 
G2-MSSM, a model which arises after considering moduli stabilization within the framework of M 
theory compactifications. A basic discussion of decay rates and branching ratios for the moduli and 
gravitinos in this model follows, with a detailed calculation left for Appendix B. Then in section 5, 
we consider again the cosmology of moduli presented in section 2) for the case of the G2-MSSM. In 
section 6, after a review of dark matter production in both the thermal and non-thermal cases, we 
consider the dark matter abundance arising from the non-thermal decay of the G2-MSSM moduli. 
This section is a more technical overview of the salient features of dark matter production, leaving 
an even more detailed treatment for Appendix A. In this section we present our main result, 
which is that the G2-MSSM naturally predicts a relic density of Wino-likc ncutralinos of about 
the right magnitude in agreement with observation. This is followed by a detailed discussion of the 
results obtained and how it depends on the qualitative (and quantitative) features of the underlying 
physics. We then conclude with considerations for the future. 

2. Ecirly Universe Cosmology in the Presence of Moduli 

Before considering the particular case of moduli in the G2-MSSM, we first briefly review the early 
universe evolution of moduli and the associated cosmological issues that can result. This section 
will also serve to set our conventions. 

Currently, the only convincing model leading to a smooth, large, and nearly isotropic Universe 
as well as providing a mechanism for generating density perturbations for structure formation is 
cosmological inflation. At present we have very little understanding of how the "inflationary era" 
might arise within the M theory framework. In what follows, therefore, wc will assume that adequate 
inflation and (p)reheating have taken place and focus on the post-reheating epoch. We will also 
conservatively take the inflationary reheat temperature to be near the unification scale 10^^ — 10^^ 
GeV, so that possibilities for high-scale baryogenesis exist. We will comment more on this issue at 
the end. 
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During inflation, the moduli fields are generically displaced from their minima by an amount 
of 0{mp) [41]. This can be seen by looking at the following generic potential experienced by the 
moduli: 

VW - Imlfti^ - %r - Hlfii, i;of + ^ ^0)^+'" (2.1) 

z rrip 

where ipo is the true vacuum-expectation-value (vev) of the field, i.e. in the present Universe. Only 
the first term in (2.1) comes from zero-temperature supersymmetry breaking, the other two highlight 
the importance of high-scale corrections and the mass-squared parameter (~ —H^^^) which results 
from the finite energy density associated with cosmological inflation [41]. As argued earlier, the 
potential (2.1) is dominated by the last two terms during inflation since Hinf ^ nisoft ^ ^3/2- 
Thus, a minimum of the potential will occur near: 

\ !/("+!) 

{i^)inf ~ -00 + mp ) H > msoft- (2.2) 

Here, for simplicity, we have implicitly assumed that the induced mass-squared parameter for 
during inflation is negative and of 0{Hf^^). This is possible for a non- minimal coupling between the 
inflationary fields and the moduli, a generic possibility within string theory. A largo displacement of 
moduli fields is also possible when the induced mass-squared parameter during inflation is positive, 
but much smaller than \Hf^j\. In this case, large dS fluctuations can drive the moduli fields to 
large values during infiation. Therefore, independent of details, the assumption we make is that 
gauge singlet scalar fields like moduli (and meson fields in the G2-MSSM) will be displaced from 
their present minimum by large values. 

After the end of inflation and subsequent cosmological evolution, when H < m-^ji^ the soft 
mass term in the potential will dominate and we have: 

{'^l') present ^ tpO H < m.off (2-3) 

V'o is also typically of order nip. In Section 4, we will present the soft masses and decay rates for the 
moduli arising from soft SUSY breaking in the G2-MSSM low-energy effective theory relevant in 
the present Universe. Thus, we see that by considering moduli in the early universe with high-scale 
inflation, it is a rather generic consequence to expect moduli to be displaced from their low-energy 
(present) minimum by an amount: 

= l(V'>m/ - Wpresentj » TOp ( ) ^ nip (2.4) 

\ '^P / 

2.1 Addressing the "Overshoot Problem" 

The evolution of moduli after the end of inflation is governed by the following equation: 

dV 

ij + {3H + r^)i;+—=0. (2.5) 
where the modulus decay rate ip XX is given by: 

r^=£'^— |, (2.6) 

which reflects the fact that the modulus is gravitationally coupled (F^ ~ Gjv ^ m~^) and is a 
model dependent constant that is typically order unity. After the end of inflation, the Universe is 
dominated by coherent oscillations of the inflaton fleld and if ~ ^ . After the decay of the inflaton 
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and subsequent reheating at temperature Tr, the Universe is radiation dominated and H ^ In 
both these phases, the evolution of the moduU can be written as: 

^ + O(l)i^+|K=0. (2.7) 

where we have neglected as it is planck suppressed. The minimum of the potential now is time- 
dependent due to the time dependence of the Hubble parameter. The evolution of the moduli in the 
presence of matter and/or radiation as in the case above, has been studied in [42-51]. In this case, 
as the modulus begins to roll down the potential, it was shown in [42,45-47]) that the presence 
of matter /radiation has a slowing effect on the evolution of the field. This can naturally allow 
for the relaxation of moduli into coherent oscillations about the time-dependent minimum^. This 
'environmental relaxation' can then slowly guide the modulus to the time-dependent minimum. 

Another possibility arises if the minimum of the potential is located at a point of enhanced 
symmetry where additional light degrees of freedom become important. This naturally arises in 
SUGRA theories that arc derived from string theories, where an underlying knowledge of the UV 
physics is known [49,50,52,53]. If the modulus initially has a large kinetic energy, as it evolves 
close to the point of enhanced symmetry, new light degrees of freedom will be produced and then 
backreact to pull the modulus back to the special point of enhanced symmetry. This simple example 
of 'moduli trapping' is present in a large number of examples in string theory with points of enhanced 
symmetry [48-51]. 

The above effects lead to a natural solution of the so-called 'overshoot problem' [54] (see also 
[55]), as argued below. As the universe expands and cools, the Hubble parameter {H) decreases 
until it eventually drops below the mass of the modulus m.^ (^ "'3/2)- Thus, from (2.1), we see 
that the first term in the potential now becomes of the same order as the other two terms and 
can no longer be neglected. At this time the modulus field becomes under-damped and begins to 
oscillate freely about the true minimum ^0 with amplitude f,^, ~ (m^ 771^)^/^"+^). As an example, 
for n = 1, fip is {irip m-ipY^'^ leading to a potential value V ^ fn^f'^ ^ m^nip which is much smaller 
than the overall height of the potential barrier at this time (~ m^m^, as in any soft susy breaking 
potential). Thus, there is no overshoot problem. 

The modulus will now quickly settle into coherent oscillations at a time roughly given by 
tosc = 2_ff^^ ^ 2m^^. After coherence is achieved, the scalar condensate will then evolve as 
pressure-less matter"^, i.e. Pm{iosc) = ^"^^/|- Because the condensate scales as pressure-less 
matter ~ its contribution relative to the background radiation 1/a* will grow with 

the cosmological expansion as a{t) ~ l/T. Thus, if enough energy is stored initially in the scalar 
condensate it will quickly grow to dominate the total energy density. 

3. Overview of Results 

This section reviews the main results of the paper without technical details. 

As explained above, the moduli start oscillating when the Hubble parameter drops below their 
respective masses. Then they eventually dominate the energy density of the Universe before de- 
caying. Within the context of G2-MSSM models, the relevant field content is that of the MSSM 
and TV -I- 1 real scalars. N of these are the moduli, Xk, of the G2-manifold and the remaining 
one is a scalar field, <j), called the meson field, which arises in the hidden sector dominating the 
supersymmetry breaking. A reasonable choice for N would be 0(50) - O(IOO). 

■^We thank Joe Conlon and Ncmanja Kaloper for discussions on this approach. 

^If there are additional terms that contribute to the potential (besides the soft mass), then a coherently oscillating 
scalar does not necessarily scale as pressure-less matter. 
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The masses are roughly as follows. The lightest particles beyond the Standard Model particles 
are the gauginos. In terms of the gravitino mass, mz/i, their masses are of order /t 7713/2, suppressed 
by a small number k. k is determined by a combination of tree level and one-loop contributions 
which turn out to be comparable. The tree-level contribution is suppressed essentially because cj) 
dominates the supersymmetry breaking, and to leading order, the gauge couplings are independent 
of ^. The precise spectrum of gaugino masses is qualitatively similar, but numerically different, 
to AMSB models. The LSP is a Wino in the G2-MSSM, similar to AMSB models. The current 
experimental limits on gauginos require that the gravitino mass is at least 10 TcV or so. In the G2 
framework, gravitinos naturally come out to be of 0(10 — 100) TeV [38]. 50 TeV is a typical mass 
that we consider in this paper. The MSSM sfermions and higgsinos have masses of order 7713/2, 
except the right handed stop which is a factor of few lighter due to RG running. Of the N moduli, 
one, Xn is much heavier than the rest, Xi. The heavy modulus mass is about 600 7713/2, while the 
(A^ — 1) light moduli are essentially degenerate with masses ~ 27773/2. Finally the meson mass is 
also about 27773/2. The decays of the moduli and meson into gravitinos will therefore be dominated 
by the heavy modulus X^. 

The decays can be parameterized by the decay width as, 

Tx=Dx^ (3.1) 
7772 

reflecting the fact that the decays are gravitationally suppressed. Dx is a constant which we 
calculate to be order one for the moduli but order 700 for (p. So, the light moduli have decay 
widths of order 10~^'^eV, corresponding to a lifetime of order 10^'^ s. The heavier scalars have 
shorter lifetimes, 10^^ s for </> and 10~^° for X^, see tables 1 and 2. So, as the Universe cools 
further and H reaches a value of order Fjc^, the heavy modulus decays. When this happens, the 
Universe is reheated to a temperature, roughly of order (rxjv"^p)^^* ~ '^^ GeV. The entropy 
is increased in this phase, by a factor of about 10^°. This greatly dilutes the thermal abundance of 
gravitinos and MSSM particles produced during reheating (by the inflaton). The abundance of the 
light moduli and meson are also diluted. Then, when H reaches order the meson decays. This 
reheats the Universe to a temperature Tj- ~ 100 MeV and increases the entropy by a factor of order 
100. Finally, as the Universe cools again and reaches a temperature of about 10""'^^ eV the light 
moduli decay. They reheat the Universe to a temperature of about 30 McV and a dilution factor of 
about 100 again. After this, all the moduli have decayed and the energy density is dominated by the 
decay of the light moduli. Since the final reheat temperature is well above that of nucleosynthesis, 
BBN can occur in the standard way. 

Furthermore, since the entropy increases by a total factor of about lO^'', the gravtino density 
produced by moduli and meson decays is sufficiently diluted to an extent that it avoids existing 
bounds from BBN from gravitino decays. 

Since the energy density is dominated by the decaying light moduli, the relic density of Wino 
LSP's is dominated by this final stage of decay. The initial density of LSP's at the time of production 
is such that the expansion rate is not large enough to prevent self-interactions of LSP's. This is 
because 

. "iff 

(3.2) 



^initial 



r^LSP > 



where the right side is to be evaluated at the final reheating temperature and av is the typical Wino 
annihilation cross-section ~ 10~''GeV~^. Therefore, the Wino's will annihilate until they reach the 
density given on the R.H.S., which is roughly 10^^ gV"^ - an energy density of 10^'^ eV*. Here we 
have assumed, as is reasonable, that since there is a lot of radiation produced at the time of decay, 
the LSPs quickly become non-relativistic by scattering with this 'background'. Since the entropy 
at the time of the last reheating s ~ lOT^ ~ 10^'^ eV^, the ratio of the energy density to entropy, 
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is around 1 eV. This should be compared to the observed value of this ratio today, which is 3.6 h"^ 
eV, where the Hubble parameter today is about 0.71. 

Therefore, we see that the Wino LSP relic density is very reasonable in these models. The rest 
of this paper is devoted to a much more precise, detailed version of this calculation. 

3.1 Scalar Decay and Reheating Temperatures 

Here we collect some more precise formulae for the decay and reheat temperatures as a function of 

the moduli/meson masses. 

The temperature at the time of decay can bo found using 

^^d - —^Sd - \ — ] g*s{Td) Ta, (3.3) 



nip V 45 



-2 _2 \ 1/3 



where = n^/s is the comoving number density and 

s='-±^ = '-^9.sT\ (3.5) 

is the entropy density with g^.^ the number of relativistic degrees for freedom^. Parameterizing the 
decay rate as above, i.e. = D^m^/ni^ we find 



For later use wc also note that if more than one modulus dominates at the time of decay then the 
temperature at the time of decay becomes 

where the sum is over all moduli (including the one that decays). When the modulus decays, the 
relativistic decay products will reheat the universe to a temperature, 

3if2 = lEl = m;^ (^^) g^T.) T,^ (3.8) 
^ Tr= (^) ' 9:'^\Tr),/f^, (3.9) 



or 



T.'[-,) ».'-W(^J ■ (310) 

Instead, if more than one modulus contributes to the energy density before decay the reheat tem- 
perature becomes 



*We will take g«s = g*, which is true if all particles track the photon temperature. This is a good approximation 
for most of the history of the universe (prior to decoupling) [56] . 
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where the sum is over all moduli (including the one that decays) and we note that this could 
lead to a subdominant radiation density compared to that of the remaining moduli. The entropy 
production is characterized by (assuming that A 1) 

where T4 and Tj. are the decay and reheat temperatures, respectively. Making use of (3.10), (3.12), 
and (3.6) we find 

A _ 2 (OKn 2\i/4 / g*s{Tr) \ ( g*a{Td) \ m^Y^ 

- - {2h0i^^f" gl'\Tr) (-^\'\^,. (3.13) 



15 \D^m^ 

For the case that more than one modulus dominates the energy density before V decays, we have 
instead 

1/4 



A= ^(2507r^) ^/%y^T.) 



/ \ 1/2 


'J2i '^^Yi' 




m^Y^ 



Y^, (3.14) 



where the sum runs over all moduli that contribute to the energy density (including the decaying 
modulus V)- 

3.1.1 Moduli decay and BBN 

Prom (3.14), we see that the decay of moduli can produce a substantial amount of entropy. There- 
fore, if any moduli present do not decay before the onset of BBN the resulting entropy production 
when decay occurs could rcsiilt in devastating phcnomcnological consequences. However, another 
possibility is provided if the late-time decay of the moduli reheat the universe to temperatures 
greater than a few MeV. Such reheating will then allow BBN to proceed as usual. Requiring that 
the modulus decay exceeds this temperature one finds from (3.10) that > 10 TeV. 



4. Summary of Results for the G2-MSSM 

In this section, we give a brief summary of the results obtained in [37-39] for the G2-MSSM. 
Readers interested in more details should consult the references above. M theory compactifications 
on singular G2 manifolds are interesting in the sense that they give rise to A/" = 1 supersymmetry 
in four dimensions with non-Abelian gauge groups and chiral fermions. The non-Abelian gauge 
fields are localized along three-dimensional submanifolds of the seven extra dimensions whereas 
chiral fermions are supported at points at which there is a conical singularity. In order to study 
phenomenology concretely one has to address the issues of moduli stabilization, supersymmetry 
breaking and generation of the Hierarchy between the Electroweak and Planck scales. These issues 
can be fairly successfully addressed within the above framework. 

In [37-39] , it was shown that all moduli can be stabilized generically in a large class of M theory 
compactifications by non-perturbative effects. This happens in the zero- flux sector, our primary 
interest, when these compactifications support (at least two) non-abelian asymptotically free gauge 
groups. Strong gauge dynamics in these non-abelian (hidden sector) gauge groups gives rise to 
the non-perturbative effects which generate a moduli potential. When at least one of the hidden 
sectors also contains charged matter, under certain assumptions defining the above framework, 
supersymmetry is spontaneously broken in a metastable de Sitter vacuum which is tuned to the 
observed value. In the minimal framework, the hidden sector, including its moduli and hidden sector 
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matter, is described by A/" = 1 supergravity with the following Kahler potential K, superpotential 
W and gauge kinetic function / at the compactification scale (~ Munif): 



N 

K/ml = -3 ln{A7r^/^Vr) + 4>(f>, ^7 = H 

i=l 

N 

/i = /2 = /hid = ^» z^; Zi = U + isi. (4.1) 

i=l 

Here V7 = ^",7^^ is the volume of the G2 manifold X in units of the eleven- dimensional Planck 
length Zii, and is a homogenous function of the Si of degree 7/3. A simple and reasonable ansatz 

therefore is V7 = YliLi s"' with positive rational numbers subject to the constraint X^iLi = 
(j) = det(QQ)i/2 = (2QQ)i/2 is the effective meson field (for one pair of massless quarks) and 
P and Q are proportional to one loop beta function coefficients of the two gauge groups which 
are completely determined by the gauge group and matter representations. The normalization 
constants C'l and C2 arc calculable, given a particular G2-manifold. /1.2 are the (trcc-lcvcl) gauge 
kinetic functions of the two hidden sectors which have been taken to be equal for simplicity, (which 
is the case when the corresponding two 3-cycles are in the same homology class). Si are the N 
geometric moduli of the G2 manifold while ti are axionic'' scalars. The integers Ni are determined 
from the topology of the three-dimensional submanifold which supports the hidden sector gauge 
groups. 

If volume of the submanifold supporting the hidden sector gauge theories (V^) is large, the 

potential can be minimized analytically order-by-ordcr in a l/V'g expansion. Physically, this ex- 
pansion can be understood as an expansion in terms of the small gauge coupling of the hidden sector 
- (ao)hid, which is self-consistent since the hidden sectors are assumed to be asymptotically free. 
The solution corresponding to a metastable minimum with spontaneously broken supersymmetry 
is given by 



where Peff = -Plii(Ci/C2). The natural values of P and Q are expected to lie between 0{1) and 
0(10). It is easy to see that a large P^s corresponds to small a for the hidden sector 

(ao i)hid = Im(/Md) « 2T,{Q-Pf ''^ ^^'"^^ 

implying that the expansion is effectively in . The (j) dependence of the potential at the minimum 
is essentially 



4 ,14 o\ 1^12 , 2 , 7 \ 



(4.5) 



Therefore, the vacuum energy vanishes if the discriminant of the above expression vanishes, i.e. if 

28(Q-P) 

= 3(g-p)-8 



^These essentially decouple from the moduli stabilization analysis. Hence they will not be considered further. 
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The above condition is satisfied when the contribution from the F-term of the meson field (F^) to 
the scalar potential cancels that from the —Sm^^^ term. In this vacuum, the i^-term of the moduli 
Fi are much smaller than F^. Since phenomenologically interesting compactifications only arise for 
Q~ P = 3 which corresponds to Peff = 84 from (4.6), we will restrict our analysis to this particular 
choice. 

4.1 Moduli Masses 

Since in this paper we arc interested in the evolution of the moduli (and meson) fields, it is important 
to study their masses in the vacuum described above. The set of gauge-singlet scalar fields includes 
N geometric moduli Sj associated with G2 manifold and a hidden sector meson field ^. Since these 
moduli and meson will mix in general, the physical moduli correspond to mass eigenstates. The 
mass matrix can be written as: 

{m%).. = (^{aiajf/^Ki+SijK2) mj/^ (4.7) 
{m').^ = {a,y/'Ksml/, (4.8) 
M0^ = ^4m^/2 (4.9) 

where Ki to K4 are obtained in [38] : 

16 \f Q Vp!s 



3 V1323y \Q-3 
32 f Q Yp^% 



Ks= \h 7^ TT^ -f^ (4-12) 



where (p^ ~ 0.734. The special structure of the mass matrix allows us to find the eigenstates 
analytically. There is one heavy eigenstate with mass tuxm = + K2Y^'^m^j2, {N — 1) 

degenerate light eigenstates with mass mxj = (i^'2)^^^^7^3/2 and an eigenstate with mass = 
{K4, — ^)^/^m3/2- These mass eigenstates of the moduli fields are given by: 



Xj = I . °'^+' — (y^<54_&1^5s' ) ; j = l,2...N-l 

Xn = xUyZv^K (4.14) 




where Ss', = \ hr^ ^^j ^-rc the canonically normalized moduli fields. The normalized moduli fields 



- V 23? ^ 

can be related to the eigenstates by (5s- = UijXj, in which Uij can be constructed using the 
eigenstates listed above. It is easy to show that {Xj)i = Uij for the eigenvector Xj. In addition, 
there is another eigenstate X^ corresponding to the meson field. Actually, the heavy eigenstate Xn 
and X^ mix with each other. This mixing hardly changes the components of the eigenstate X^ and 
Xff, since tuxn ^ "^0- However, the mass of the eigenstate X^ is affected by the mixing. The masses 
itixn and only have a mild dependence on Q (for PeS = 84, Q — P = 3), and do not depend 
on the number of moduli N at all. The mass of the light moduli mx^ does not even depend on Q. 
Taking the expression for K2, one immediately finds that mXj ~ 1-96 1713/2, j = 1,- ■ ■ ,N —1. This 
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result is very important since light moduli are then not allowed to decay into gravitinos, essentially 
eliminating the moduli induced gravitino problem. Choosing a reasonable value of Q to be of 0(10), 
one finds that rricf, is roughly around 2 77x3/2 while mx^ is roughly around 600 Changing values 

of Q by 0(1) hardly changes the moduli masses mx^ and m^. Therefore, the above typical values 
will be used henceforth in our analysis. To summarize, the meson and moduli masses in the G2 
-MSSM can be robustly determined in terms of 7713/2. 

4.2 Couplings and Decay Widths 

Understanding the evolution of the moduli also requires a knowledge of the couplings of the modiili 
(meson) fields to the visible sector gauge and matter fields. Since all the moduli are stabilized 
explicitly in terms of the microscopic constants of the framework, all couplings of the moduli and 
meson fields to the MSSM matter and gauge fields can in principle be explicitly computed. Here 
we focus on the moduli couplings to MSSM matter and gauge fields. A different visible sector, as 
might arise from an explicit construction, will give rise to different couplings of the moduli fields in 
general, although with roughly the same moduli masses. 

Here we will give a brief account of the important couplings of the moduli meson to visible 
gauge and matter fields and set the notation. Details are provided in Appendix B. The most 
important couplings of the moduli and meson fields involve two-body decays of the moduli and 
meson to gauge bosons, gauginos, squarks and slepton, quarks and leptons, higgses and higgsinos. 
The three-body decays are significantly more suppressed and will not be considered. 

Let us start with the decay to gauge bosons and gauginos. The relevant part of the Lagrangian 
is given by: 

■^gauge boson, gaugino — QX^gg -^k ^ ^i/^ ftX^gg -^k A A -f- 

964,ogg -^^0 A^A^ k = 1 ■ ■ ■ N; a = 1,2,3 (4.15) 

Here, Xk, Si^o, F^v ^^id A" are the normalized moduli, meson, gauge field strength and gaugino 

fields respectively. The expression for the couplings will be provided in Appendix B. It is important 
to note that the meson field does not couple to gauge bosons since the gauge kinetic function 
does not depend on ^o- The normalized moduli eigcnstates Xk have already been discussed. The 
others can be written as: 

where fsm is the gauge kinetic function for the visible SM gauge group. In the rest of the paper, 
we will neglect the hats for these normalized fields and rup in the couplings for convenience. 

The coupling of the moduli and meson fields to the MSSM non-higgs scalars (ie sfermions) turn 
out to be important, as will be seen later. Since the Standard model fermion masses (including 
that of the top) are much smaller than that of the moduli, the decay of the moduli and meson to 
these fermions will not be considered. The coupling to the MSSM sfermions can be written as: 



.'non-higgs scalars 



'Xff 



■0// 



where fa are the canonically normalized scalar components of the visible chiral fields Ca, i.e. 
fa = ?° . The couplings to the higgs and higgsinos arc different due to the presence of the higgs 

bilinear Z HuH^ + h.c in the Kahler potential [39], which gives rise to contributions to the ji and 



- 11 - 



B/j. parameters. In addition to the couplings similar as those in Eq.(B.17), there are additional 
couplings for scalar higgses, which can be schematically written as: 



■^higgs 3 gxHdH^ XjHdHu + g'xHdH-a d^iXjd'^{HdHu) + c.c 

d^6<t>od>'{HdH^)+c.c (4.18) 

As explained in [39], all higgs scalars except the SM-like higgs and all higgsinos are heavier than 
the gravitino, implying that the moduli and meson fields can only decay in this sector to the light 
SM-like higgs {h). The coupling to the SM-like higgs can be determined from the above coupling 
as explained in appendix B. 

Finally, the moduli and meson fields can also decay directly to the gravitino. In fact, it turns 
out that the (non-thermal) production of gravitinos from direct decays dominates the thermal 
production of gravitinos in the early plasma. Therefore, it is important to consider the moduli and 
meson couplings to the gravitinos. Since the meson and light moduli are lighter than twice the 
gravitino mass (as seen from the previous subsection), only the heavy modulus can decay to the 
gravitino. 

The explicit form of these couplings in terms of the microscopic constants is provided in ap- 
pendix B. An important point to note is that these couplings are computed from the theory at a 
high scale, presumably the unification scale. However, since the temperature at which the moduli 
decay is much smaller than the unification scale, one has to RG evolve these couplings to scales 
at which these moduli decay (around their masses). The RG evolution has also been discussed in 
appendix B for the important couplings. Once the effective couplings of these moduli and meson 
are determined, one can compute the decay widths, as shown below. 

For the G2-MSSM model, we have found that light moduli and meson dominantly decay to light 
higgses and squarks, while the heavy modulus dominantly decay to light higgses only. In appendix 
B, we have explicitly calculated the decay widths of the moduli Xk and meson. The widths of 
moduli can be schematically written as: 

^^r' = 1,2^" ~ ^Xk^gg + ^Xk^gg + ^Xk^qq + ^Xk^hh (4.19) 



where k = 1 • • • iV and Nq = 12 is the number of gauge bosons or gauginos. Note that A^'' is 
significant only for fc = 1, 2...,N— 1 (see appendix B). For the meson, the width can be written as: 



= ^(^oAt + At + At) ( ^ j . (4.20) 



4.3 Nature of the LSP 



Before moving on to discuss the evolution of moduli in the G2-MSSM, it is important to comment on 
the nature of LSP in this framework. As explained in detail in [39], the G2-MSSM framework gives 
rise to Wino LSPs for choices of microscopic constants consistent with precision gauge unification. 
Therefore, in our analysis we focus on the Wino LSP case. As we will see, a Wino LSP turns out 
to be crucial in obtaining our final result. 
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5. Evolution of Moduli in the Gs-MSSM 

In this section, we apply the general discussion in Section 2 to the model of the G2-MSSM reviewed 
in the previous section. For clarity we will summarize our main results focusing on the more salient 
aspects of the physics, leaving the more technical details of the calculations to Appendix A. We 
will illustrate our computations with benchmark values, in order to get concrete numerical results, 
and comment on the choice of the benchmark values in section 7. 

As discussed in Section 2, we assume that cosmological inflation and (p) reheating have provided 
adequate initial conditions for the post-inflationary universe. 

5.1 Moduli Oscillations 

As reviewed in the last section, we have a heavy modulus Xtv, A^ — 1 light moduli X^, and the scalar 
meson </>. These will begin to oscillate in the radiation dominated universe once the temperature 
cools and the expansion rate becomes comparable to their masses. 

For a benchmark gravitino mass vahic^ of 50 TeV, the heavy modulus will begin oscillations 
first, at around t^^ « 10"^^ seconds, corresponding to a temperature of roughly T = 10^^ GeV. 
Following the heavy modulus, the other moduli will begin coherent oscillations around IQ-^" s 
corresponding to a temperature of roughly 10^^ GeV. These results are summarized in Table 1 
below. 



Modulus 


Mass (m3/2 = 50 TeV) 


Oscillation Time (seconds) 


Xn 


mx^ = 600 TO3/2 


if^^ = 2 X 10-32 


x^ 


< 2 TO3/2 


ttc = 7 X 10-30 


Xi 


"fnxi ^ 2 77x3/2 


t^,^ = 7 X 10-30 



Table 1: Oscillation times for the G2-MSSM moduli 

Since coherently oscillating moduli (pm) scale relative to radiation as Pm/Pr ~ cL{t) ~ 1/r, 
the moduli will quickly come to dominate the energy density of the universe, which is then matter 
dominated. Following the beginning of coherent oscillations of the heavy modulus, until the decay 
of all the moduli the universe will remain matter dominated. We will see that this, along with 
the entropy produced during moduli decays, results in negligible primordial thermal abundances of 
(s)particles compared with the non-thermal abundances coming from direct decays of the moduli. 
This will be crucial in addressing the gravitino problem and establishing a Wino-like LSP as a 
viable dark matter candidate through its non-thermal production. 

5.2 Moduli Decays and Gravitino Production 

As the universe continues to cool the expansion rate will eventually decrease enough so that the 
moduli are able to decay. This occurs when H ~ Tx , at which time the moduli will decay reheating 
the universe and producing substantial entropy. We will parameterize the decay rates of the G2- 
MSSM moduH as: 

rx=Dx^, (5.1) 

where Fx is the decay width for particle X. The decay times will be computed for a set of bench- 
mark values of Dx for the various moduli (meson) which can be obtained by choosing particular 
(reasonable) sets of values of the microscopic constants (see appendix B for details). 

®We give detailed numerical values for 7713/2 = 50 TeV. It will be clear that values a factor of two or so smaller 
or larger than this will not change any conclusions in this and related analyses. 
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Modulus 


Decay constant 


Decay Time (seconds) 






Tx„ = 9 X 10-11 




= 710 


= 6 X 10"^ 


Xi 


Dx, = 4.00 


Tx, = 10 X 10-^ 



Table 2: Decay constants and lifetimes for the G2-MSSM moduli for a set of benchmark microscopic values 

5.2.1 Heavy Modulus Decay and Initial Thermal Abundances 

Given the G2-MSSM values in Table 2 above, the heavy modulus will be the first to decay at 
around lO^^is. This decay will produce a large amount of entropy Axjv — Safter/ Sbefore ~ 10^° 
(even though the energy density of the heavy modulus is less than that of the meson and moduli), 
reheating the universe to a temperature T^" = 41 GeV. The entropy production will not only 
dilute the thermal abundances of all (s)particles, but also all the other moduli. One particularly 
important non-relativistic decay product of the heavy modulus is the gravitino. Gravitinos will be 
non-thermally produced by the modulus decay with a branching ratio B^'^ = 0.07%, which yields 

a comoving abundance ^3^"' ~ ~ 10^*^. This can be compared to the thermal abundance 

of gravitinos, which before modulus decay is Yy^™'^^ = 2.67 x 10"^. This is further diluted by 
entropy production resulting from the decay, i.e. Y^^^^™"''' —^ YlJ^™°-^ I bi.Xfi ~ 10"^^. We see 
that the thermal contribution to the gravitino abundance is negligible compared to that from non- 
thermal production. A similar result follows for all other (s)particlcs that arc thermally populated 
following inflation. Therefore, the primary source of (s)particles, and in particular gravitinos and 
Lightest SUSY Particles (LSPs), will result from non-thermal production resulting from decays of 
the moduli. 

5.2.2 Meson/Light Moduli Decays and the Gravitino Problem 

The decay of the heavy modulus is followed by the decay of the meson, at around 10~^s (for 
benchmark values). The meson will decay before the light moduli because of a larger decay width 
compared to that for the light moduli (see appendix B for details). Similar to the heavy modulus, 
the meson contribution to the energy density is small compared to that of the iV — 1 light moduli. 
Nevertheless, it produces some entropy (A,^ « 121) and reheats the universe to a temperature of 
around 134 MeV. The entropy production will again dilute the abundance of light moduli, and any 
(s)particlcs present, including the gravitinos from the heavy modulus decay. 

The decay of the meson to gravitinos is particularly important, as this can result in the well- 
known gravitino problem. If the scalar decay yields a large number of gravitinos, these gravitinos 
can later decay producing a substantial amount of entropy that could spoil the successes of BBN. 

The entropy produced from the decay of the meson and the other light moduli further dilutes 
the gravitino abundance from the heavy modulus. The primary contribution to the gravitino relic 
abundance comes from the decay of the heavy modulus since the other fields have masses of order 
2 7713/2. After the decay of the meson, the energy density of the N —1 light moduli is the dominant 
contribution to the total energy density of the Universe. 

Given that the A'' — 1 light moduli are approximately degenerate in mass, their decays will occur 
at nearly the same time, after the decay of the meson. The resulting reheat temperature is found 
to be approximately 32 MeV, which is an acceptable temperature for consistency with the bound 
of 1 MeV set by BBN [57-60]. 

We note that the moduli decay rates have a strong dependence on the gravitino mass (as it sets 
the moduli mass scale;). So, the decay of the light moduli being able to avoid BBN constraints is a 
result of the fact that the gravitino mass is relatively large (rn3/2 ^ 50 TeV). However, as explained 
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in detail in [39], the gauginos are significantly suppressed relative to the gravitinos allowing us to 
still obtain a light (< TeV) spectrum which can be seen at the LHC. The decay of each modulus will 
contribute to the total entropy production, and one finds that the total entropy production for the 
set of benchmark values of the microscopic constants is given by Ax^ = 418. We also note that the 
light moduli lifetime depends inversely on the decay constant , so if instead of taking relatively 
large values Dx^ = 4 we take relatively small values = 0.4, we find a reheat temperature of 10 
MeV which is still compatible with BBN^. The decay of light moduli to gravitinos is kinematically 
suppressed for the same reason as for the meson. The final gravitino abundance is then just the 
contribution from the heavy modulus decay diluted by the decay of the meson and light moduli and 
is Y^J™'' = Y^^^/Axi ~ 10~^^ . The above gravitino abundance is well within the upper bound 
on the gravitino abundance set by BBN constraints, as it will not lead to any significant entropy 
production at the time the gravitinos decay. Thus, we find that there is no gravitino problem, 
in the G2-MSSM. In addition to the relativistic decay products, the light moduli will also decay 
appreciably into neutralinos (LSPs), which we consider in detail in the next section. 

6. Dark Matter from the Ga-MSSM 

Natural models of electroweak symmetry breaking (EWSB) require additional symmetries and 
particles beyond those of the Standard Model. The additional particles typically come charged 
under additional discrete symmetries suppressing their decay to Standard Model particles (e.g. R- 
parity, KK-parity, etc.), so such models predict an additional, stable, weakly interacting particle 
with an electroweak scale mass, i.e. they naturally predict a candidate for Weakly Interacting 
Massive Particle (WIMP) cold dark matter. In the case of the G2-MSSM, this gives rise to a 
Wino-like neutralino which is the lightest supersymmetric particle (LSP) of the theory. 

For completeness in section 6.1 we will review the standard calculation for computing the 
(thermal) dark matter relic density today. In section 6.2, we will then revisit this calculation for 
non-thermal production of LSPs resulting from scalar decay. In Section 6.3, we examine how non- 
thermal production is naturally realized in the G2-MSSM and predicts the Wino LSP as a viable 
WIMP candidate. 

6.1 Standard Thermal Dark Matter 

In the standard calculation of the relic abundance of LSPs it is asstuned that prior to BBN the 
universe is radiation dominated. In particular, it is assumed that the dark matter particles are 

created from a thermal bath of radiation created from (p)reheating after inflation. In this radiation 
dominated universe, the Friedmann equation reads = m~^pr, with pr = (7r^/30)(?*T^ the 
radiation density and g'* the number of relativistic degrees of freedom at temperature T. 
The evolution of LSPs are given by the Boltmann equation 

fix = -3Hnx - (cTv) [n% - nig] , (6.1) 

where (av) is the thermally averaged cross-section, nx is the number density, and rieq is the number 
density of the species in chemical equilibrium, i.e. XX ^ 77, where 7 is a relativistic particle such 
as the photon. 

Assuming that initially the dark matter particles are relativistic {mx < T) and in chemical 
equilibrium, then they will pass through three phases as the universe expands and cools. Initially 
their density will be determined by all the factors on the right side of (6.1). As long as the 

interactions of the particles take place on smaller time scales than the cosmic expansion then the 
particles will remain close to their equilibrium distributions. While the species is relativistic [rnx < 

'^See appendix A for a discussion of the range of the coeiRcients Dx^ ■ 
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T) this means that their comoving abundance is given by Yx = nx/s « = const.. Once the 
universe cools enough from the cosmological expansion so that X becomes non-relativistic (T < mx) 
then particle creation becomes more difficult (Boltzmann suppressed) and the comoving abundance 
tracks that of a non-relativistic species Yx ~ Y^ = 0.145 x^/^ exp(— .t) where x = m,x/T. The 
particle density will continue to decrease until the number of particles becomes so scarce that the 
expansion rate exceeds the annhiliation rate and the particle species undergoes 'freeze-out'. Prom 
(6.1) we see that at this time the number density is given by: 



tlx 



{av) 



(6.2) 



eg 



where Tf indicates that this relation only holds at the time of freeze-out. Using (6.2) and Yx ~ Y-^^ 
at the time of freeze-out, we find that freeze-out is only logarithmically sensitive to the parameters 
of the model, = « \n\mxmp{av)\ and corrections are O (In In a;/). Taking both the cross- 
section and mass mx to be weak scale at around 100 GeV we find that = 4 and thus the 
freeze-out temperature is T/ = mx/25 w 4 GeV. Prom (6.2) and (3.5), we find the comoving 
density at freeze-out: 

^/ = (6-3) 



45 1 

45 1 / nix 



(6.4) 



27rVl0ao5i 



1/2 



X/, (6.5) 



where we have taken (cru) = a^m^ . We note that this answer is rather insensitive to the details 

of frcczc-out, and the abimdance is determined solely in terms of the properties of the produced 
dark matter (mass and cross-section). In particular, there is no dependence on the underlying 
microscopic physics of the theory. 

6.2 Non-thermal Production from Scalar Decay 

We know from the successes of BBN that at the time the primordial light elements were formed the 
universe was radiation dominated at a temperature greater than around an MeV. However, perhaps 
surprisingly, there is no evidence for a radiation dominated universe prior to BBN. In particular, we 
have seen that in the presence of additional symmetries and flat directions, scalar moduli can easily 
dominate the energy density of the universe and then later decay. The presence of these decaying 
scalars can alter the standard cold dark matter picture of the last section in significant ways. 

To understand this, consider the decay of an oscillating scalar condensate which decays at 
a rate ~ m^/nip. When the expansion rate becomes of order the scalar decay rate {H ~ F) 
the scalars will decay into LSPs along with relativistic (s)particles which reheat the universe to a 
temperature T^. If this reheat temperature is below that of the thermal freeze-out temperature of 
the particles Tf ^ mx/25 then the LSPs will never reach chemical equilibrium. As an example, if we 
consider a scalar mass ~ 10 — 100 TcV this gives rise to a reheat temperature ~ ^^T^nip > 
MeV where F^ m^^/rn^. The decay of </) in a supersymmetric setup could lead to LSPs with 
weak-scale masses mx ~ 100 GeV, which have a thermal freeze-out temperature Tf ~ mx/25 
few GeV. We see that in this case Tr < Tf is quite natural and the particles are non-thermally 
produced at a temperature below standard thermal freeze-out. Thus, the particles will be unable 
to reach chemical equilibrium. 

Depending on the yield of dark matter particles from scalar decay, there are two possible 
outcomes of the non-thermally produced particles. 
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6.2.1 Case one: LSP Yield Above the Fixed Point 

If the production of LSPs coming from scalar decay is large enough, then some rapid annihilation is 
possible at the time of their production. Since the particles are produced at the time of reheating, 
we know from the Boltzmann equation (6.1) that the critical density for annihilations to take place 
is: 

3H 



^ {av) 

which is different from the result (6.2) in that here the reheat temperature and not the thermal 

freeze-out is the important quantity. This is very important bccaiisc Tr ^ \/^4>m,p depends on 
the microscopic parameters of the theory as the reheat temperature is set by the decay rate of 
the scalar. In the standard case, we saw that the freeze-out temperature, or more precisely, the 
parameter Xf = mx/Tf was only logarithmically sensitive to the parameters of the dark matter and 
gave no information at all about the underlying theory from which the dark matter was produced 
(e.g. scalars from the underlying microscopic physics). 

Given that the initial number density of particles exceeds the above bound {nx{0) > n^), the 
LSPs will quickly annihilate until they reach the density (6.6). Thus, the critical value n'j^ serves 
as a fixed point for the number density, since any production above this limit will always result in 
the same yield of particles given by rix- From this one finds the comoving density [40] 

ci 1 ^^std ( 



(6.6) 

Tr 



gl/^ m,p{(7V)Tr \Tr J 

where Ci = 45/ (27r\/T0). Wc sec that non-thermal production can yield a greater comoving density 
than standard thermal production by a factor (T//Tr). For the example considered above, namely 
~ 10 — 100 TeV, mx ^ 100 GeV, and Tr ~ few MeV we find the comoving density is enhanced 
by a factor ^ 10^ — 10"^. One interesting consequence of this is it allows room for larger annihilation 
cross-sections for the LSPs. For example, in standard thermal production a Wino-like LSP leads to 
too small a relic density since its annihilation cross section is only s-wave suppressed . In the case 
of the G2-MSSM, non-thermal production is a natural consequence of the microscopic physics and 
a Wino LSP will provide a perfectly suitable WIMP candidate. 

6.2.2 Case two: LSP Yield Below the Fixed Point 

The other possibility is that the decay of the scalar yields few enough LSPs (nx(0) < n^^) so that 
annihilation does not occur. Then the comoving abundance is simply given by 

>0. = i^,A-rr~^, (6.8) 

where is the branching ration of scalars to LSPs and Y^^^ is the initial abundance of scalars in 
the decaying condensate. We note that again this result depends on the underlying physics of the 
UV theory, since both the branching ratio and the reheat temperature are coming from the physics 
of the scalar. 

6.3 Dark Matter in the G2-MSSM 

As shown in [39], the LSP in the G2-MSSM is predominantly Wino-like. There are two significant 
sources of these LSPs in the G2-MSSM - direct production from decays of both the gravitino and 
the light moduli. As explained earlier, the thermal abundance of LSPs in the early plasma after 
infiation is vastly diliitc;(l by the entropy productions from the heavy modulus, meson and the light 
moduli. Therefore, the thermal abundance of LSPs is negligible. In addition, the LSPs produced 
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from decays of the heavy moduhis and the meson field are also diluted by the entropy production 
from the light moduli and are negUgible as well. 

The light moduli may decay to LSPs directly, or via decay to superpartners. Prom Section 

4 the branching ratio for this process to occur for a set of benchmark values of the microscopic 
paramaters is B^^p ~ 25% and the comoving abundance is then found to be: 

y^sp = ^xMhp^N - l)4f = 1.19 X 10-^ (6.9) 

where = 417.7 [(A'')/100]^^^ is the entropy production from the decay of all the light moduli Xi. 
Here we have taken benchmark value for the number of light moduli to be 100. The corresponding 
number density at the time of reheating is 

riLSP = s{T^')Ylsp. (6.10) 
= 1.79 X 10"" GeV^ (6.11) 

As discussed in the last section, we must compare this number density of LSPs to that of the critical 
density for annihilations (6.6). At the time of reheating from the light moduli the Hubble parameter 
is given by 

^ 5* \ \-^r ) _ A AO in-22 



^(^^^n^J ^=4.48xlO-GeV. (6.12) 
The dominant (s-wave) annihilation cross section for the LSPs (W°W^° W~^W~) is given by 



(av) = aomfsp = -J- = 3.26 x lO"^ GeV-^ (6.13) 

''''LSP " \ ~ 

where Xw = i^%,/fn\spi = 80.4 GeV is the W^-boson mass, and g2 ~ 0.65 is the gauge coupling 
constant of SU{2)l at temperatures Tr ^ McV, and this defines ctq. It is crucial that the cross- 
section is s-wave so that there is no temperature dependence in {(Jv). We will comment more on 
this in section 7. Using (6.6) we find the fixed point density for annihilations 

nlsp = 4.12 X IQ-^^ GcV^. (6.14) 

We sec that the produced density is greater than the fixed point value n'^gp > n^gp and annihila- 
tions will occur. This corresponds to the "LSP yield above the fixed point" case discussed above. 
Thus, the LSPs produced will quickly annihilate down toward the fixed point value in less than a 
Hubble time. The relic density of dark matter is then given by the fixed point value (6.14) and the 
critical density of LSPs today coming from decay of the light moduli is 

^x ^ m.spYEsp ^ 1 ( 45 \ f m^^A ^ ^ 

Pc/so Pel SO \2TTy/m§:ao) \mpT^^^ ) ^ ' 

where sq and pc are the entropy density and critical density today, respectively, and we have used 
the experimental value pc/so = 3.6 x l{)~^h? GeV with h parameterizing the Hubble parameter 
today with median value h = 0.71. 

In addition to this contribution, there is also the contribution from the decay of non-thermal 
gravitinos produced from the heavy modulus which have a final abundance Y^J^"''' w 10""^^. The 
contribution from gravitinos to the critical density of dark matter is then 

, 7- final 

^(3/2) ^ 612. ^ 0.0008/l-^ (6.16) 

Pc 

which is negligible compared with that coming from the light moduli. 
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Thus, the total critical density in dark matter coming from the LSPs of the G2-MSSM is : 

o ^2 n^-/ misp \3 / 10.75 y/-* / 3.26 X 10-^GeV-2\ / 4 V amgp V 100 TcvV^'' 

"^^^'^ ~°-'UMGivU^J I W Jl^J l^^J ' 

(6.17) 

where we have included all the parametric dependence of the answer derived in Appendix B. This 
value should be compared to the experimental value VLcDMh? = 0.111 ± 0.006 [61]. For those used 
to {av) in other units, note that 1 GeV~^ = 0.4 x 10~^^cm^. 

This result is not presented in terms of central values - rather it is the best value we can obtain. 
The LSP mass can be larger than 100 GeV, but not smaller. The decay constant Dx^ can be order 
4, but a scan of the microscopic parameter space suggests a somewhat smaller value for the only 
calculable example so far known (see appendix B.4). A better understanding of the string theory 
could give 4 or a larger value. Whereas 7713/2 is somewhat constrained to be at most about 100 
TeV by the parameters of the framework, as explained in [39]. Therefore, this framework is rather 
constrained and predictive. We view the closeness of this result as a success, and as an indication 
that improving the underlying theory may improve the agreement with data. 

7. Discussion of Results 

We have seen in the previous sections that for natural values of microscopic parameters, there is no 
moduli and gravitino probem in realistic G2 compactifications. In addition, within the G2-MSSM, 
the non-thermal production of Wino LSPs from the light moduli give rise to a relic density with 
the right order of magnitude (up to factors of a few). It is possible that with a more sophisticated 
understanding of the theory, one could obtain a result more consistent with the observational 
results. It is also worthwhile to understand these results from a physical point of view. The results 
obtained above depend surprisingly little on many of the details of the microscopic parameters. 
In particular, there is essentially no dependence of the final relic density on the total number of 
moduli (A'^), the masses {mxNt^4>) a^^d couplings {Dx,^,Drj,) of the heavy modulus and meson 
fields as well as the initial amplitudes of the moduli (fxk) and meson (/^) fields. This is good 
in a sense since our understanding of the underlying theory and many of the above microscopic 
parameters is incomplete. However, the result does depend crucially on certain qualitative (and 
also some quantitative) features of the underlying physics, as we discuss below. In general it is 
better if results depend on the microscopic theory, since then data can tell us about the underlying 
theory. 

One very important feature which helps avoid the gravitino problem is that the meson and light 
moduli have masses which are of order (actually slightly below) two gravitino masses, as we saw 
explicitly in Section 4.1, This kinematically suppresses their decays to the gravitino. The gravitino 
abundace is thus dominated by decay of the heavy modulus which is further diluted by entropy 
production from the decays of the meson and light moduli. Therefore, a natural mechanism for 
solving the gravitino problem in a generic setup is that the modulus which decays last does not 
decay to the gravitino. The moduli problem can also be easily solved in frameworks where the 
gravitino mass is > 10 TeV, which is naturally satisfied in the G2 framework. 

Another qualitative feature of the G2 framework is that there is a hierarchy in the time scales 
of decay of the various moduli (meson) fields. Since the mass of the heavy modulus is much larger 
(^ 300 times) than that of the other moduli (meson), it decays much earlier. Also, from our current 
understanding of the Kahler potential of the meson and moduli fields, one finds (see appendix B) 
that the meson decays before the light moduli due to a larger decay width. The precise computation 
of the decay width depends on the nature of the Kahler potential for the meson and moduli and the 
Kahler metric for matter fields, and one might argue that there are inherent uncertainties in our 
understanding of these quantities. However, the only qualitative feature relevant for cosmological 
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evolution is that the meson decays before the hght moduh. As long as the light moduli decay last 
(which we have argued in the appendix to be the natural case from our current understanding of 
the Kahler potential), the result does not depend on any of the masses and couplings of the heavy 

modulus or the meson field. The final result depends only on the masses and couplings of the light 
moduli which decay last. The same qualitative feature could be present in other frameworks arising 
from other limits of string/M theory. 

Now that it is clear that it is the light moduli decaying at the end which affect the final relic 
density, it is important to understand their effect more closely. In any theory of (soft) supersymme- 
try breaking, the mass of the light moduli will be set by the gravitino mass scale. In the context of 
low energy supersymmetry, the gravitino mass will typically be in the range 1 — 100 TeV. Therefore, 
the light moduli will also be typically in the above range^. Since the reheat temperature of the 
moduli basically depends on the moduli masses (assuming the coefficient is C(l)), the light 
moduli will typically give rise to a reheat temperature T^^ of 0(1 — 100) MeV, which is far smaller 
than the freezeout temperature of the LSPs (T^^^ ~ GeV) which could be produced from the light 
moduli. This is true for the G2 framework and could be true for many other frameworks as well. 
Therefore, with T^* < T^^^ , the final outcome for the relic density will depend on the whether 

the number density of the LSPs produced from the light moduli {n\^gp) is greater or smaller than 
the critical number density at T^' ("^isplT-^i )• 

For the G2 framework, for natural values of the microscopic parameters one finds that n^^p > 
n^gplyXj as shown in section 6.3. This is equivalent to the inequality: 

^LSP^Xi > — ( o ) « 1^07 
(To ^X- 

with 7 = (7.1) 

'^3/2 

where (Tq is defined by (6.13) and we have used mx^ ~ 1.96 7713/2- As explained in [39], the quantity 
7 depends predominantly on 5, which characterizes the threshold correction to the gauge couplings 
at the unification scale. The dependence on other microscopic parameters such as V7 and C2 (see 
section 4) is largely absorbed into the gravitino mass. The suppression factor 7 depends almost 
linearly on \5\, and typically takes value in the range ~ (1 — 6) x 10~^. Now, the constraint (7.1) 
is easy to understand. For natural values of microscopic parameters in the G2 framework, one has 
B^^gp = 0{25%), Dxi = 0{1) (see appendix B) which easily satisfy (7.1) above. In order for other 
frameworks to realize this situation, a criterion similar to (7.1) needs to be satisfied. 

When (7.1) is satisfied, the final relic density can be written as (see (6.14) and (6.17)): 

2 _ rnLSpY£sp _ ^ ( 45 \ / mlgp 



^Lsph ~ ^^^^^ ^^^^^ l^27rV10g,croy' [m^^T^ 

1 / 45 



. U GeV-'^ = 1« ""V-"" {^\ ,7.2) 

An upper bound on the observed value of the relic density implies that smaller values of 7 and 
niLSP and larger values of Dx^ are preferred. A small 7 implies that for a given LSP mass a 



*This is however not true for Large Volume compactifications as the lightest modulus in that case is much lighter 
than [25]. 
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heavier gravitino is preferred implying that the moduh be correspondingly heavier. Also, since 7 
is roughly linear in \6\, smaller values of \d\ are preferred. These features can be seen easily from 
the plots in figures 1 and 2. Figure 1 shows a contour plot of the relic density in the Dx^ — m3/2 
plane for two (large and small) values of \S\ which correspond to two (large and small) values of 7. 
Figure 2 shows the dependence of the relic density on the reheat temperature of the light moduli 




25 50 75 100 125 150 175 200 



m3/2 (TeV) 

Figure 1: The contour plot of the relic density in the G2-MSSM in the Dxi — ^13/2 plane for two (large 
and small) values of \S\ which correspond to two (large and small) values of 7. The solid lines are for 5 = —3 
(a correction to a"^^^ of order 3/26), and the dashed Unes for S = —4.5. 

(r/^'). As seen from the first line in (7.2), the relic density is inversely proportional T^^ implying 
that a higher reheat temperature is preferred. A higher T^' corresponds precisely to a larger Dx^ 
and (larger 7713/2) explained above. 

As explained in section 6.3, the nature of the LSP is also crucial to the final result for the relic 
density. For the G2 framework, the annihilation cross-section is s-wave and does not depend on 
T^' . On the other hand, if the LSP were Bino, the cross-section would be p-wave suppressed and 
would depend linearly on T^' /misp- thereby making it suppressed relative to the s-wave result. 
This would make the relic density much larger than the result obtained for the s-wave case above. 
Therefore, the upper bound on relic density prefers small mixing angles (or vanishing mixing angles, 
as in the G2-MSSM) with the Bino and Higgsino components. This can be seen from figure 3. 

8. Summary and Future Directions 

In this paper we have emphasized the importance of the cosmological moduli and gravitino problems 
and the relation to adequate generation of dark matter in thermal equilibrium, or generation of too 
much dark matter non-thermally in string/M theory frameworks. Focussing on G2 compactifica- 
tions, in particular on the G'2-MSSM, we have found that the decay of moduli in this framework is 
rather naturally consistent with BBN constraints, and the associated large entropy production at 
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Figure 2: The LSP relic density for the G2-MSSM plotted as a function of the reheat temperature of 
the light moduli. The solid lino assumes no coannihilation with charged Winos; the dashed line includes 
coannihilation with charged Winos. 

late times (but before BBN) results in an avoidance of the gravitino probleni(s). Moreover, we have 
seen that the late decay of the light moduli into Wino-like neutralinos leads to a nearly acceptable 
relic density of cold dark matter. This result arises from a combination of entropy production and 
LSPs from moduli decay giving an adequate relic density from non-thermal production of dark 
matter. This process offers an explicit example of how thermal dark matter production is not the 
dominant source of cosmological dark matter, especially in the presence of moduli. The LSP is 
Wino-like here as well as in anomaly mediated theories, but for interestingly different reasons - 
here the tree level gaugino masses are universal but about the same size as the anomaly mediated 
ones, and the finite one loop Higgsino is comparable with both. 

The result for the final relic density depends parametrically on the couplings and mass of the 
light moduli (which decay last) and the mass of the LSP. The masses of the light moduli and the 
LSP are set by the gravitino mass scale and depend on a set of underlying microscopic parameters 
of the theory. The couplings of the moduli depend on the Kahler potential of the theory. Since 
our understanding of the Kahler potential is incomplete, it is only possible to make reasonable 
assumptions to proceed, which is what we have done, but one can see that most of the results are 
insensitive to these uncertainties. That is because the moduli decays produce a large number density 
of LSPs, which then annihilate down to the final relic density that only depends on the reheating 
temperature. Prom (7.2) and figure 1, we see that an upper bound on the relic density prefers a light 
LSP, heavy gravitino and large couplings to the visible sector parameterized by Dxi (defined in 
(4.19)). These results obtained have been explained in terms of the underlying qualitative features 
of the framework. These qualitative features could be present in other string/M theory frameworks 
as well, leading to similar results. 

There is not yet a satisfactory inflation mechanism for the G2-MSSM. This is under study. 
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Figure 3: The LSP relic density in the G2-framework plotted as a function of the mixing angle of Bino 
and Wino for M2 = lOOGeV. 

Fortunately, our results are not sensitive to that. We assume only that at an early time inflation 
ends and the energy density of the universe is dominated by moduli settling into the minimum of 
the potential. 

In future, one would like to understand the origin of the baryon asymmetry in the Universe 
(BAU) within string/M theory frameworks. In the G2 framework, the large entropy production 
resulting from the decay of the moduli was crucial for addressing the gravitino problem. However, 
this entropy will also act to reduce any initial baryon asymmetry. Therefore, one requires a large 
initial asymmetry or a late-time mechanism for regeneration of the asymmetry. For example, a large 
initial baryon asymmetry could arise from the Affleck-Dine mechanism [62] , or it could happen that 
the superpartner parameter space allows for late-time electroweak baryogenesis. This is work in 
progress. 

Understanding the above issues would be crucial to solving the "cosmological inverse problem" 
(see [63,64] for some preliminary work in this direction), usually considered separate from the "LHC 
Inverse Problem" [65]. Within the context of realistic string/M theory frameworks, however, the 
two inverse problems merge into one "inverse problem" as the microscopic parameters characterizing 
the underlying physics of any framework have predictions (at least in principle) for both particle 
physics as well as cosmological observables, thereby providing unique connected insights into these 
basic issues. 
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A. Cosmology of the G'2-MSSM Moduli - A detailed treatment 

In this appendix, we include detailed calculations leading to the abundances, entropy production, 
and reheat temperatures quoted in the paper for sets of benchmark values of the microscopic 
parameters. The computation of couplings and decay widths of the moduli and meson fields in 
terms of the microscopic parameters which motivate the benchmark values will be given in appendix 

B. We have retained the parametric sensitivity to the gravitino mass, number of moduli (topology), 
and the overall couplings of the moduli (meson) in order to address the robustness and plausibility 
of the framework. 



A.l Heavy modulus oscillations 

At the time the heavy moduli {X^) starts coherent oscillations the universe is radiation dominated 
and the Hubble equation is given by 

3H^ = 3( ^mxS = fj^") g.T\ (A.l) 



The temperature at which the modulus starts oscillating is then given by 

^osc \ j \ osc )i 

228.75 \ ^'^ ( mx, ^ 



= 2.70 X 10^2 """•''^ GeV. (A.2) 



Prom this we find the entropy density 



27r^ 

= 1.98 X 10^9 f ( ^ GeV^, (A.4) 

V 228.75/ UoOm.voy ^ ' 



and the comoving abundance is then 



==^'"^./i»»-'(Ti"). 



4.51 X lOM iSAV'Vii^V f25i=HiV'\ (A.5) 



(Tosc) J J \ nix. 

The oscillating modulus will quickly come to dominate the radiation density and the tempera- 
ture at this time is given by 

= 1.80 X 10- ( { { ' GeV, (A.6) 



\ 9* {Tosc )) V 600 m3/2 ) V "ip / 

so that we see once the modulus starts coherent oscillations it quickly overtakes the energy density 

rXt. 
eq 
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A. 2 Meson and Light Moduli Oscillations 

Because the meson and light moduli are approximately degenerate in mass (i.e. = rrijc.) they 
will begin to oscillate at the same time, 

= 3 (\n,^ = m-^ {^^9*{TL) {tL^ + mx^Yx. s{TtS) ■ (A.7) 

Noting that the radiation term has already become negligible compared to the heavy modulus 
density we find the temperature at this time is given by 



2 2 ll/3 

4> P 



\ 1/4 / \ 2/3 \ 1/6 / \ 2/3 

228.75 \ I \ /SOOma/aN I rup ^ 



(A.8) 



= 8.24x101° [^^\ —^'^ ^ GeV (A.9) 

which is in excellent agreement with the exact answer obtained numerically (including radiation) 
^otc = 9-97 X 10^°. The entropy density at this time is 

... . (^) " [^y (^) ^ o.v^ 

The meson (p initial abundance is then 

^(0) _ . ,n . 1 n6 r 228J5 \ ^^V /, \ V fx. \ ' f mx. \ '''(2 m,,. 



y]"^ = 5.30 X loM [^\ ^ — ^ • (a.h) 

\9*{Tosc)J \mpj \mpj \600m3/2j \ J ^ ' 

The light moduli will begin coherent oscillations at roughly the same time as the meson. Their 
abundance is then given by 



y(°) = (iv-i)yW 



N-l\ f 228.75 y-'^ ffxA' ffx.^Y ( ^x. \'''' (2ra^i^ 



5.25x10 )\g*{To,c)) V mp J \ ) VeOOmg/aj V mx. J' ^^'"^^^ 

where we have implicitly assumed that because the masses of the meson and light moduli are 
approximately degenerate they will have equal oscillation amplitudes^. 

A. 3 Heavy Modulus Decay 

Once the Hubble parameter decreases to the point when H w ^x. , the heavy modulus decays and 
from (3.11) the corresponding reheat temperature is, 

=41.40 (^Y' (B^Y' f^^Y' (^Y' fmY' GeV. (A.13) 

To understand the N and 4> dependence in this expression, we note that from (3.11) the reheat 
temperature includes the factor, 

1/4 



Using that the meson and light moduli have degenerate mass and therefore equal oscillation ampli- 
tudes (i.e. mXiYxi = {N — l)m^Y^) we find 



1+N ^ « iV ^ (A.15) 

mx.Yx.J V rnx.Yx./ 



^We note that initially this may not be the case, but at the onset of coherent oscillations (much less than a Hubble 
time) the system will settle into this symmetric configuration. 
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which leads to the parametric dependence in the reheat temperature. 

Using (3.14) the entropy increase resulting from the heavy modulus decay is 



Axjv = 4.35 X 10 



io^54r.^")\'^V 228.75 \"'(fxA" ( 1-6 



10.75 ; \g,{T^,P:)J \rnpj \Dx 



X 



mxff J \mpj \100 



(A.16) 



where we have again used (A. 15). Therefore, after the decay the other moduli abundances are given 
by 



1.22 X 10- i^^T i^X''^ i^Il^Y' (A.17) 



^ 1.21 X 10- i^Y' (l^Y" (^!^] ' (A.18) 

where we have again used N — 1 ^ N. There is also a decay to gravitinos with branching ratio 
^3/2^ = 0.2% = 0.002. The corresponding comoving abundance is thus, 

y(0) 

^3/2 - ^ ^ -"3/2 X -T— , 

= 1.45 X 10- (^Y' c^Y^ f^ii^Y'^ (^Y^ ^^)^'^(A.19) 

\0.07%J \g,(T^'^)J \ 1.6 J \600m3/2j \U) V JV 7 

A. 4 Meson Decay 

When the meson decays, its contribution to the total energy density will be less than that of the 
other N — 1 light moduli. The universe will be matter dominated before and after the decay, but 
because the two energy sources are comparable there is a somewhat significant entropy production. 
The meson decay reheats the universe to a temperature 

\N J \g4Tr)J \711.6j \2ms/2j ^ ' 

The entropy increase is given by 

* V 1-6 / \ ) V600m3/2y \ J \mpj ^ ' 

The decay of the meson will further dilute the other moduli, we find 



— i4)"'(^)"'(^r(^)'" 

The decay of both the meson and the light moduli to gravitinos is kinematically suppressed, so that 
the only source of gravitinos comes from the decay of the heavy modulus. This abundance after 
the decay of the meson is then 

-''"3/2 ~ '^<j> -'"»3/2 ' 

1.19 X 10-" I ^ i'^Y^ (^\''' (J^Y^ i^Y^ 

0.07% \N J \g.{Tr)J \ni.6j \2m3/2 J \ mx^ )\U) 
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A. 5 Light Moduli Decays 

The decay of the light moduh results in a reheating temperature 

which agrees with the bounds set by BBN (i.e. T^* > 1 MeV). The resulting entropy production is 



The new gravitino abundance is given by 

2.86 X 10-^^ I ^] ('-^) f^Y' (^Y'^ C-!^] (R^Y' 

0.07% J \ N J \g4Tr)J \2ms/2j \ mx^ J \ ^ J \ f4, J 

which is small enough to avoid the gravitino problem. The light moduli will decay into LSPs 
yielding an abundance 

^LSP ~ '-^Xi^LSP^ Xi ' 



where B^^p is the branching ratio for the decay of the light moduli to LSPs. This corresponds to 
a number density at the time of decay of tilsp = 1-79 x 10~^^ GeV^. 

As wo noted in the text, this abundance is produced below the freeze-out temperature of the 
LSPs (non-thermal production) and is greater than the critical density (6.6) for annihilations to 
take place, which is rix. = 4.12 x 10"^^ GeV'^. Thus, the LSPs will quickly annihilate (in less than 
a Hubble time) and the final abundance will be given by the critical value. 

Thus, the relic density coming from the decay of the light moduli is given by 

Pc 

= 0.26/^- f^I^Y (lO^Y^ (^.2..1U -'GeV-^X /^^/^ ^jn^yf^ / lOOWX^^^ 
VlOOGeV/ \9*{Tr)J \ <Jo ) \DxJ \ mx, J \ mg/a / 

where sq and pc arc the entropy density and critical density today respectively, and we have used 
the experimental value Pc/so = 3.6 x lQ~'^h? GeV where h parametrizes the Hubble parameter 
today with median value h = 0.71. 



B. Couplings and Decay Widths of the Moduli and Meson Fields 

In this section, we discuss the moduli couplings to MSSM particles and then calculate their decay 
widths in terms of the microscopic parameters of the G2-MSSM framework. This will motivate the 
benchmark values used for numerical results throughout the paper. We will find that the moduli 
decay into scalars is very important. 

B.l Moduli Couplings 

Let us first consider the couplings associated with A'' eigenstates Xj of the geometric moduli s^. For 
simplicity, we neglect the small mixing with the meson modulus cj) (we will return to that later). 
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First consider the moduli coupling to gauge bosons through the gauge kinetic function The 
relevant term is: 

£ D -ilm(/,„)F;,F''''- (B.l) 

= - l(Im(/,„))F;,F«''- - 1 ^ NrSsiF^^^F'^^'^ (B.2) 

i 

where we have expanded the moduli as Si = {si) + Ssi. After normalizing the gauge fields and the 
moduli fields, the interaction term can be written as: 



AT 



where B and Cj are defined as: 



i;^E^rJ^c/., x,f;^f'^^-^ (b.3) 



6V2 



i=l 

V7 



N \ -1 



\i=l 



^ E ^ (B.6) 



For the heavy modulus, since (-^jv)f = fcii, we have Cn = |B ^ while for the light moduli 
Xi, i = 1, - ■ ■ , (A'' — 1), it is easy to show: 



^C! = fsmH, (B.7) 

i=l 

where I is the length of the vector Xj^ defined as {Xj^)i = {Xiq)i Nf"^/Ni and is the angle 
between X^ and Xn- So, generically Cj are less than one. There are two extreme cases: one when 
N?^ = kNi in which the moduli couplings to gauge bosons vanish since the vector Xn is orthogonal 
to Xj, and the other when Xp^ equal to one of the Xi's in which all Cj's are zero except one. 

For the couplings to gauginos, the dominant contribution comes from the following terms in 
the lagrangian: 

C D —difsruF'X'X' + h.c. (B.8) 

where difsm = AT?'" and —i arises because of the convention of the moduli chiral fields Zi =ti + ist 
we used. Expanding the F-terms of the moduli fields around their vevs, we have: 

F' = {F') + {ds,F')6sk (B.9) 

The derivative of the F-term can be calculated as follows: 

= -le -m3/2 [^J^N.u ^ + --(3a, - 25,,)u^ + —{3a, - 25,,) j 

= -ie-'-'ms/a [-^SiNkhib2V - S^a, + 26,, + • • • ^ (B.IO) 
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where in the last hne, the subleading terms are not exphcitly shown. 7 is the phase in the super- 
potential which will be set to zero for simplicity without affecting any result here. We have used 
the following equations: 

- - 2 - 

After normalizing the moduli fields and the gauge fields, the couplings are given by: 

AT , AT 



(B.ll) 



1 



XiXX 



1^3/2 



fc=l •'^'^ fe=l ^^'^ 



XiX^X" + h.c. 



14 



12 



-^3/2 



XiX^X" + h.c. 



For the light moduli fields, the first term vanishes and the couplings turn out to be: 



■'XiXX 



B Cim3/2 XX + h.c. 



For the heavy modulus field, the first dot product is unity and the coupling is: 



XiXX 



12 



^3/2 (^^'^1^2) XiA"A° + h.c. 



(B.12) 



(B.13) 



(B.14) 



The moduli couplings to other MSSM particles can be derived generically by expanding all the 
moduli around their vevs in the supergravity lagrangian: 



C D ks.pVj*'^V^~f + iK^pp^a^Vj" - V{rJ) + ■■■ 



(B.15) 



where /" and /" are fermions and their superpartners. The other derivative terms involving moduli 
and matter fields are not explicitly shown for simplicity. The relevent coupling here are the moduli- 
sfermion-sfermion coupling and the moduli-fermion-fermion coupling. They are found to be 

C D ds,K^0 \Ssi dj*^d^~f + iSsi P^a^d^f] - d,,m% 5si~n~f + ■■■ 



la 

- 9 Xiff 



d^{Xi jnd^t + c.c. + iX, f^a'^dj^ - g'^j^.XJrf^ + 



(B.16) 
(B.17) 



where /" and /" arc the canonical normalized fields. For simplicity, we consider the Kahler metric 
to be diagonal Ka/s = Kaiap, then 



la 

9 xjf 



2s^ 



14 

^"^3/2 



(B.18) 

(B.19) 



where {X'l^)f = £,i,a{^N)i/ cii and ^i^a = Sidg. \og{Ka). In this calculation, we have used the fact 
that d^gKa = and have neglected terms involving F-terms of geometric moduli which are 
suppressed relative to 7713/2- 

For the coTipIings to the higgs doublets, there are differences from other scalars. The kinetic 
terms and the mass terms for the higgs fields in the MSSM can be written as: 



C D Kh. 



df,H*d''Hu + iHua''df,Hu\ +■■■ 
- {Bn'HdHu + C.C.) 



(B.20) 



-29- 



where 



(B.21) 



is only generated by the higgs bihnear term in the Kahler potential [39]. To derive the modular 
couplings to higgs doublets, one needs dg-fj,' , which is: 



(B.22) 



One can see that the second and the third terms are of order 7723/2 while the rest are suppressed. 
Therefore, the dominant contribution is: 



1 / 4 

dsilJ'' ~ ^'rn3/2ida^Z) i- -SmNibib2i' + 4:Si, 



(B.23) 



For simplicity, taking all the phases of the superpotential and that of Z to be vanishing, we find: 



(B.24) 



gXjHuHu 



m 



'3/2 



+ log Kh, 




2 72 
™3/2-^eff 



Xisr ■ Xj + 4X^' • Xj 



[X'l,)"^ . XA + ^mi,^{X'l,r- ■ X, 



14 



(B.25) 



where {X'^)i = ^(Xjv)i and Q = Sidg. log(Z). We also use the fact that d^^^Z = and the i^-terms 
Fi/mp ^ 7773/2 for geometric moduli. To get the corresponding couplings for H^, we can simply 
replace iJ„ by in the above equations. The coupling of moduli to higgs through the kinetic term 
is similar to the non-higgs scalar 



la 

9 XiH^Hu 



{X'^)"" ■ Xi 



(B.26) 



Let us now consider the B/j, term, which is given by: 

B^J.' = (2771^/2 + Vo)Z - my^F'^drnZ + m3/2F"'[d„,Z - dm log{K„^KH,) Z] 
- F^F^[drndnZ - dn\0g{kH^kH,)drnZ]. 



The corresponding derivative is given by: 
1 



dsiBfi' « ^mli^Zds. \og{KH„KHd) [-■^SiNkbib2i' + 25,^ 



'2tJI^ l2^Si Z ^ 



which gives rise to the coupling: 

-C D gxjH^H^XjHdHu + ex. 
\/l4 , 



-^™3/2 



Z^sl^'l^'b^b^i^f^^C-^^N-Xj 



2(X;^)^" • X, + {Hu ^ Hd) + 4X™ • X, 



(B.27) 



(B.28) 



(B.29) 



(B.30) 
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Besides the term mentioned above tliere is another coupHng from the bihnear term in the kahler 
potential K ^ Z{si)HdHu + h.c. [40]. This term leads to a coupling: 

^ 3 g'x.H,H^d^Xjd^{HdHu) + ex. (B.31) 
= —^Z^sX'l^-Xj (B.32) 

This coupling could be very important since it is proportional to the moduli mass squared if equa- 
tions of motion of Xi are used. Again for the coupling to be unsuppressed, the bilinear coefficient 
Z should have a sizable dependence on the geometric moduli Sj, which is natural. This coupling is 
essential for electroweak symmetry breaking in the G2-MSSM. 



B.2 Meson Couplings 

In the G'2-MSSM framework, the hidden sector is not sequestered from the visible sector and there 
are couplings between the hidden sector meson field and various MSSM particles, which we want 
to compute. First since the tree level gauge kinetic function does not depend on there is no 
coupling to gauge bosons. However there are couplings to the gauginos which depend on d^„, which 
are computed to be 

4s- 

d^,F' = -ie-'^— ^.Fm3/2, (B.33) 
300 

^=^+2+|+0(P,i^). (B.34) 

After normalization of fields, the coupling of meson to the gauginos is given by: 

jO-s^oW = e"*'^— i— •?^TO3/2^?!'oAA (B.35) 

We now move on to the couplings of the meson field to scalars. We will assume that the Kahler 
metric and the higgs bilinear Z do not depend on 0o- We then have for the non-higgs scalars: 

1 Otti ^ ~ ~ 
= V2m3/2 (500 '>nz/2)5(t>oI* I 

« V2ml,^M^ + ^)5M*f (B.36) 

In the above, we have neglected terms proportional to Fi/nip which are <C TO3/2- There are various 
kinds of couplings of the meson to the Higgs fields and H^. The coupling originating from 
the term J d^O {ZHuHd + c.c) does not give rise to any contribution since Z is assumed to be 
independent of The couplings jCs^igH^H^ and CsipoH'Ha are computed as follows: 

1 d{k^%T+m'l) 



V2{Z^f, + l)m' 



3/2V0 



(B.37) 
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^5(t>oH*Hd ^^"^ ^® obtained from the above by replacing i7„ with Hd- Again, wc have neglected 
terms proportional to Fi/nip. Finally, we look at the coupling Cs^^h^h^- It is given by: 



1 d{Bfi') 



N 



(B.38) 



The coupling Cs^qHih* can be computed by taking the complex conjugate of the above expression. 



B.3 RG evolution of the couplings 

In the last subsection, we computed all the relevant couplings of the moduli and meson at a high 
scale, presumably around the unification scale. However, since the scale at which moduli decay 
is much smaller than the unification scale, one should in principle use the efii'ective couplings at 
that scale to compute the decay widths. The RG running of the moduli-scalar-scalar couplings are 
especailly important for the third generation squarks and the higgs doublets and are the main focus 
of this subsection. The leading contribution to the /3 functions are terms proportional to and 
5f3^°, which are given below: 



1 



X0Q3Q3 



3|2/t| gXjHaH^ 
|2„/ 



0{9XjH^H„) « Y^^'^*'^ {gXjH^H^ +Xt) , 

Pig'x.H^Hj ~ Y^^\yt\^9'x,H^H^^ 

9x,Q,Qs (2|2/tl' - ySl) +2\yt\'Xt 
1 , 16 



16 



167r2 

1 

167r2 



gx.usus [ 8\yt\' - y53 ) +^\ytrXt 



gx. 



I |2 16 2 

ml - yff3 



(B.39) 



where Xt = gx^H^H^ + 9x Q3Q3 + fl'^j«3«3- For other beta functions not listed above, the RGE 
eff'ects can be neglected. 

To examine the RG effects on the moduli-scalar-scalar couplings, wc take all the weighted dot 
products involved in the moduli-scalar-scalar couplings to be equal for simplicity^ ^, 



X'n ■ Xi — {X'^)" ■ Xi — II. 



(B.40) 



^'^Here we have not included the digrams proportional to gx^gg and QXjgg, since their contributions are relatively 
smaller 

^^The more general case will be studied later. 
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This is reasonable as their structure is very similar. So the high scale couplings can be written as: 



to,H„H„ = 9X,H,H, = ^™3/2(3^eff + (B.41) 

6 



g'x.H^H^ = 9x,H,H, = (B.42) 



4V14 



gxjH^H^ = — g— m^/2^eflFn (B.43) 



\/14 

9X,H,H^ = V^effH (B.44) 

Using the beta functions given in Eq.(B.39), wc can see that at low scale gXjH^H^ is squashed 
because of the large yukawa couplings. Similarly 9xjQsQ3 ^^'^ 9XjU3U3 decrease significantly and 
become negative at low scales. 

One important thing to compute for moduli decay to light higgs is the effective coupling 9'^xjhh-' 
which can be written in terms of the couplings to higgs doublets 

dxfhh = {9XiH„H„ - '^mlg'x.H^H J cos'^ a + {gx.H^H^ - "^rnlg'^.H^Hj sin^ a 

- [gXiHaH^ - mx.g'x.H^H J sin 2a (B.45) 

where all the couplings involved should be evaluated at low scales and a is the higgs mixing angle. 

For the G2-MSSM, the higgs sector is almost in the "decoupling region", which implies a « /3 — 
■|. Now with universal boundary condition for the weighted dot products for concreteness and 
simplicity, the effective coupling of moduli to hh final state is given by: 



9xih ~ ^"^3/2 [(3^eff + l)(sin' a + Ki cos^ a) - 2i^2^eff sin(2a)] U (B.46) 

where Ki and K2 are the RG factors. To estimate these factors, we take yt = 1, ct^nif ~ ^^-^ 
Zes = 1.58, which is the same as the first Benchmark G2-MSSM. Then, typically we find Ki ^ 0.2 
and K2 ^ 0.5. For readers not familiar with the details of the G'2-MSSM, it is helpful to know 
that generically tan/? ~ 1 and Z^s ~ 1-5. For the effective coupling to third generation squarks, 
including the RG effects, we have: 



4i.3-<.ig3-4^-3/2n (B.47) 
where q!/^^ = Qy if — nt^q' From the above RGE results, we find that the couplings to 

^Xjff ^X-iff f^Xjff ' ^ 

the non-higgs scalars and higgs should be roughly of the same order because of the large radiative 

correction even when some of them arc suppressed relative to the other at the high scale boundary. 
Therefore, if the couplings to scalars are large, then we should expect a significant branching ratio 
of the moduli to LSPs. 

For the coupling of the meson field to scalars, the /3 functions are exactly the same. Similar 
to the analysis of light moduli, wc introduce factors Ki and K2 to account for the RG effects on 
g<t>H^H^ and g,pH^H„- Typically one has Ki ~ 0.25 and K2 ~ 0.5. From Eq.(B.37) and (B.36), we 
find the coupling g^H^H^ is at least Z'^^J^ ~ 30 times larger than g^^j at the high scale. Because 
of this large coupling g^H^H^i even if the couplings g^q^q^ and g,f,z^ui are zero at the high scale, 
they can still be generated at the low scale, which is proportional to g^iH^H^ by a factor K3 ~ 0.1. 

B.4 Decay Rates of the Moduli 

Now that we have computed all the the relevant couplings for moduli decay, we can proceed to 
compute the corresponding decay widths. In the following, we give the result of decay widths for 
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all the moduli, calculated from the two-body width formulae. There could be contribution from 
three-body decays, which is generally small because of the phase space. Although certain three- 
body decays, e.g. moduli to top quarks and higgs [8,40] is relatively large, it is still comparatively 
small in the current framework compared to the two-body decay modes. 
For light moduli Xi,i = 1, ■ ■ ■ , {N — 1), the total decay width is 



Dx rn\r 7 / v v v v \ ni^- 

T{X,) = — = ^ (NgA^^ + NgA^^ + Af^ +A^^)^, (B.48) 
nip I ZTT V / rrip 



where Ai*, A^*, A^* and A^' are the corresponding coefficients for the decays to gauge bosons 
gg, gauginos gg, non-higgs scalars // and light higgs bosons hh respectively. They are given by: 

(B.49) 
(B.50) 
(B.51) 




1 2 



(3Z,^ff + l)(sin^a + i^i cos a) — 2K2Zef[ sin(2a) 



(B.52) 



Here, weighted dot products in the scalar couplings are assumed to be equal and are denoted as 11 
as in the last subsection. In addition, the RGE effects on the couplings are included. In the above 
result, the gaugino and gauge bosons are treated as massless. The two-body decay to the standard 
model fermions is suppressed by i^^)'^ ^ 10^^, so it is neglected in our result; even the top 
quark contribution is small. For the decay to non-higgs scalars, naively there is a large kinematic 
suppression since these scalars have mass close to m3/2. However, the RGE running significantly 
decreases the third generation squark mass at the scale much lower than the unification scale. In 
G2 MSSM framework, the lightest stop is tn which is about 4 times lighter than the gravitino. It, 
therefore, has a large contribution to the partial width. In addition, Q3 {bL and ii) are also light 
enough such that they contribute to the decay width. 

The above result for ^'s depend on the specific choices of the fundamental parameters, such as 
Ui, Ni and iV/™, through several weighted dot products of vectors Xn and Xi. These quantities 
are different for different moduli. However from Eq.(B.7) they are constrained by: 

JV 

J2{X'n ■ Xif = IX^vl' sin^ e. (B.53) 

i=l 

Similar constraints apply for other products. From the above equation, one expects that on average 

{X'^-X,f-j^\X'^?sm^e (B.54) 

which is suppressed by 1/{N — 1). It is obvious that this symmetric configuration is favored in 
cosmology. If one wants the moduli to decay before BBN, then the most dangerous modulus is the 
one with the smallest total decay width, which is bounded by the average width. This gives rise 
to a strong constraint on the geometry of the G2 manifold, since the width is suppressed by the 
number of moduli N. In the following discussion, we will focus on this symmetric configuration. 
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In order to evaluate the decay width and the branching ratio, one needs to know the typical 
values of these weighted dot products of Xn and Xi. To do the estimation, we generate a set of 

and Q randomly with the following conditions: 



fundamental parameters a.i, Ni, iVf™ 



N 

1=1 



7 

3' 



1 < iVf" < 2, 2<N,<6, -1< 6 < 0, -1< < 0. 



(B.55) 



The above ranges are chosen based on constraints arising from the G2 framework and our current 
understanding of the Kahler metric of visible matter fields in realistic constructions. We also 
impose the supergravity condition V7 > 1 and volume of three-cycle Vg™ w 26. The ranges of Ni 
and N^"^ are chosen such that the efficiency of the parameter generation is maximized when the 
above constraints are imposed. Due to our primitive understanding about the kahler metric, the 
modular weights (corresponding to and Q) are taken randomly in the allowed range. We plot 
the distribution for B~^{X'j^ ■ Xi)^ and (Xj^' • Xi)^ in Fig. 4, where we can see the typical values are 
2 X 10-4 and 20. This resuh can be understood from the very rough estimate (X^ ■ Xi) ^ ^Jai ~ 



l/y/N and {X'^ ■ X,) - l/7a7 ^ VN. The distribution of {{X'j^)°' ■ X.f is expected to be about 
the same as [X'^ ■ Xif ^ since they all have the same structure. However, one should be aware that 
all the weighted dot products are independent and so are not necessarily equal. 




0.0001 0.0002 0.0003 0.0004 0.0005 0.0006 



B"2 Cl^ 




(Xn .Xi)' 



Figure 4: Left: distribution of the average of B''^{X'i^ ■ Xif. Right: distribution of the average of the 
weighted dot product {X'^ ■ Xi)^. 



Now let us estimate the decay width for the light moduli. Consider the first benchmark model 
of G2-MSSM [39] for example, assuming the weighed dot products take their average value, we find 
~ ^ IQ-^, ^3 ' ~ 7.3 and ^4 ' ~ 20.5. To summarize, the main channels of interest for 
light moduli decays and their partial widths are T{gg) = T{gg) « 0.024 sec"^, r(t/jt/{) « 60 sec""'^, 
r(?L?L) = r(6L^L) ~ 43 sec~^ and T{hh) w 412 sec^^. The total width is the sum of these partial 
width. LSPs arise mainly from gauginos (including LSPs), tt and bb, so the LSP branching ratio 
is the sum of the gaugino and squark channels divided by the total width. One can see that the 
decay to higgs and scalar dominate the decay of the light moduli. The total decay width is about 
558 sec^^ or the corresponding = 0.86. The branching ratio of the light moduli to LSP is 
about 26%. These results should still be roughly correct for other benchmarks, differing at most 
by 0(1) since the dependence on the mass spectrum is mild as seen from the explicit result of Af^. 
The main uncertainty arises from the deviation of those weighted dot products from their typical 
values. To explore the more general case, one can relax the condition that all the weighted dot 
products are equal. Instead we choose: 

X'^ ■ X, = {X';,)"- ■ X, - {X'i,)''- ■ X, = Hi (B.56) 
{X';,)'^' ■ X, = {X';,f' ■X,^U2 (B.57) 
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Then we vary Hi and 112 according to the distributions of the weighted dot products in Fig. 4. The 
distribution for Dxi and the branching ratio to LSPs is shown in Fig. 5. One can see that the 
branching ratio has a very small variation, but the distribution of Dxi has a long tail. In the paper, 
we will use 0.4 < < 4 for concreteness, although other values may be possible. 



[ 




. 3 



Figure 5: Left: distribution of Dx^- Right: distribution of moduli branching ratio to LSP. 



For the heavy modulus X^, the total decay width is 

T{Xn) = — ( A^G^f " + iVc^f " + " (B.58) 

( ZTT \ / TOp 

where ^f"" corresponds to the decay to gauge bosons gg, gauginos gg and higgs bosons, and are 
given by 

= 4 (B.59) 



9 \mx^ 

Zlsix'^-X^"" ■ (B.61) 



In the above result, we have not included the contributions from the decay to non-higgs scalars and 
fermions since they are suppressed by {rn^/2/mXfiY ^"^^ ("^//'^^w)^ given the large mass of the 
heavy modulus mx^ ~ 600 x 7713/2- Taking benchmark 1 of G2-MSSM in [39] and typical values for 
weighted dot products, we get ylf" w 0.1, " w 0.01 and " « 50. The total width is about 
3 X 10^° sec~^ or the corresponding Dx^ ~ The branching ratio to LSPs is about 3 x 10"'^. 

The decays of moduli to gravitinos is also very important. The decay of a modulus to gravitinos 
can be calculated using the following formula: [19] 

|^(e//)|2 5 

where q'"^^^'^ is the effective coupling of the modulus field to gravitinos which includes effects of 
moduli mixing. For the heavy modulus, the coupling arises from the mixing with meson field, 
since the goldstino is mainly the fcrmionic partner of the meson. Since the heavy modulus is much 
heavier than the meson a rough estimate^^ gives Gx'^'^'^ ^ "^3/2/"^X]v Therefore, for the heavy 



^^There could be an additional suppression in special cases as discussed in [7, 12]. We thank Fuminobu Takahashi 
for discussions regarding this point. 
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modulus, the decay rate to gravitino is 



T{Xn ^ 2^-3/2) 



28877 ml 



(B.63) 



This corresponds to B. 



3/2 



7 X 10- 



In addition, since the heavy modulus decays much earlier 



than other moduli, the gravitino produced will be diluted by the subsequent moduli decays. So 
this estimate is enough for our discussion of gravitino problem. For both the light moduli and the 
meson fields, the decay to gravitino is kinematically suppressed since mxi,m^g « 2m3/2. 

B.5 Decay Width of the Meson 

The total decay width of meson modulus is: 



r(^.^o) = = ^ {NgA^ + At + At 



mi, 



(B.64) 



where Af corresponds to the decay to gauginos gg, non-higgs scalar // and light higgs bosons hh, 
and are given by: 



At 



A. 



00 _ ( 

2(f>l V ^0 

a \ 

1/2 



2, 2^^ V 




N 



2 2J- 



sin 2a 



(B.65) 



(B.66) 



(sin^ a + Ki cos^ a) 



(B.67) 



Here, as discussed in the last subsection, the low scale couplings to third-generation squarks are 
dominantly generated from RG running and are related to the coupling g^iH^H^ by a factor ~ 0.1. 
For the first benchmark of G2-MSSM in [39] and taking the simplest assumption = Qi = —1/2, we 
get Af" « 13.7, At « 2.9 x 10^ and At « 1-3 x 10^. One can see that this result is enhanced from 

the naive estimate by the total number of moduli ^ TV and large (hidden-sector) three-cycle volume 
v. The total decay width is about 4.6 x 10^ sec^^, corresponding to = 711. The branching ratio 
to LSPs is about 18%. Again, this can change by 0(1) for other benchmarks. 
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